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ABSTRACT 



Context. This paper is the second in a series devoted to studying the properties of binaries with M giant primaries. 
Aims. The binary frequency of field M giants is derived and compared with the binary fraction of K giants. 

Methods. Diagrams of the CORAVEL spectroscopic parameter 5*6 (measuring the average line width) vs. radial-velocity standard 
deviation for our samples were used to define appropriate binarity criteria. These then served to extract the binarity fraction among 
the M giants. Comparison is made to earlier data on K giant binarity frequency. The Sb parameter is discussed in relation to global 
stellar parameters, and the Sb vs. stellar radius relation is used to identify fast rotators. 

Results. We find that the spectroscopic binary detection rate among field M giants, in a sample with few velocity measurements 
(~ 2), unbiased toward earlier known binaries, is 6.3%. This is less than half of the analogous rate for field K giants. This difference 
originates in the greater difficulty of finding binaries among M giants because of their smaller orbital velocity amplitudes and larger 
intrinsic jitter and in the different distributions of K and M giants in the eccentricity - period diagram. A higher detection rate 
was obtained in a smaller M giant sample with more radial velocity measurements per object: 11.1% confirmed plus 2.7% possible 
binaries. The CORAVEL spectroscopic parameter Sb was found to correlate better with the stellar radius than with either luminosity or 
effective temperature separately. Two outliers of the Sb vs. stellar radius relation, HD 190658 and HD 219654, have been recognised 
as fast rotators. The rotation is companion-induced, as both objects turn out to be spectroscopic binaries. 

Key words, binaries: spectroscopic - stars: late-type 



1. Introduction 

This paper is the second in our series discussing the 
spectroscopic-binary content of a sample of M giants drawn 
from the Hipparcos Catalogue (ESA 1997), for which 
CORAV EL radial v elocities ha ve been obtained in a system- 
atic way dUdrv et al.|[T997: Famaev et alJlIOOS,) . The main driver 
behind such an extensive monitoring effort lies of course with 
the kinematics; indeed, the kinematical properties of the present 
samp le of M giants have been fully analysed by iFamaev et alj 
(l2005l) . But this large amout of material may also be used to 
search for b inaries, which is disc ussed in a series of three papers. 
In Paper I dFamaev et alj|200"9l) . we present the radial-velocity 
data and orbital elements of newly-discovered spectroscopic bi- 
naries, and also discuss the intrinsic variations sometimes caus- 
ing pseudo-orbital variations. As a follow-up, the present paper 
presents the first attempt to derive the frequency of spectroscopic 
binaries in such an extensive sample of M giants. A side topic of 
this paper is the identification of fast rotators among M giants, 
using the CORAVEL line-width parameter Sb. This way of find- 



* Based on observations carried out at the Swiss telescope installed at 
the Observatoire de Haute Provence (OHP, France), and at the 1.93-m 
OHP telescope 

** Postdoctoral Researcher, F.N.R.S., Belgium. Currently at 
Department of Physics, Technion-Israel Institute of Technology, 
Haifa 32000, Israel 

Postdoctoral Researcher, F.N.R.S., Belgium 
Research Associate, F.N.R.S., Belgium 



t 



ing binaries is independent of radial velocities. Finally, we inves- 
tigate the relation between the CORAVEL line-width parameter 
Sh and stellar parameters and conclu de that Sh is best cor related 
with the stellar radius. In Paper III (iJorissen et al.ll2009l) . more 
evolutionary considerations are presented, based on the analysis 
of the eccentricity - period diagram of our sample of M giants. 

2. Binary frequency among M giants 

The sample used for the present study and its three subsam- 
ples are extensively described in Paper I. Sample I consists of 
77 1 stars from the Hipparcos survey, for which at least 2 radial- 
velocity measurements have been obtained (Udrvetal. 1997|). 
Sample II consists of about one-third sample I stars (254 stars), 
for which at least 4 measurements have been obtained. Sample 
III consists of 35 sample II stars with cr(K-) > 1 km s~^, and 
sample IV of 138 sample II stars with (T{Vr) < 1 km 

2. 1. Spectroscopic binaries (sample I) 

As explained in Paper I, the detection of spectroscopic binaries 
(SBs) among M giants cannot rely solely on a test compar- 
ing the radial-velocity standard deviation cr(V^r) to the average 
instrumental error < e >, because of the mass motion existing in 
the atmospheres of these stars (either due to convection or to pul- 
sation), causing intrinsic radial-velocity variations. In Paper I, 
it was indicated that the parameter Sb, measuring the average 
width of spectral lines (connected for the instrumental width). 
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provides a useful tool to distinguish intrinsic radial- velocity jit- 
ter from orbital motion (see Fig. [TO This is because the Sb (i.e., 
line-width) parameter is a good tracer of the evolutionary state of 
giant stars (as is shown in Sect. 12. 3t . and extensive mass motions 
are occurring in the envelopes of highly evolved giants. These 
mass motions cause the radial velocity jitter, which is thus large 
in highly evolved giants, where Sb is the largest. 

In the Sb — ctq diagram, where ao{Vr) — (ct^— <e>^)^/^, 
stars are flagged as binaries if they fall above a 'dividing line' 
(Fig ID- The binaries identified from the literature and from the 
more intensive monitoring of subsamples II, III, IV (Sect. 3.1 
and Table 4 of Paper I) have been used to define that bound- 
ary: the vast majority of these binaries fall above that boundary, 
while at the same time, confirmed non-binary stars are not found 
in that region. A dividing line with a slope of 0.57 fulfills these 
conditions. There is one apparent exception in the lower left cor- 
ner, HD 40282, which in follow-up observations turned out not 
to exhibit orbital motion. It may be a statistical fluctuation in 
the measured jitter (Compare to the discussion below, where a 
Gaussian distribution for the radial velocity jitter is assumed.) 
Besides, its status as an exception depends on the exact defini- 
tions of Sb and of the dividing line. It must be stressed that this 
slope for sample I is steeper than that obtained for the dividing 
line in diagrams where the radial-velocity standard deviation is 
computed from more extended data sets characterising samples 
II, III, and IV (see the corresponding panels in Fig.!!). 

As these lines are defined by the observed distribution of the 
single-star objects, the difference in the position of the dividing 
lines can be explained by a combination of two factors: the num- 
ber of radial-velocity data points for a given object and the total 
number of objects in a sample. 

First, every ao value is an observational estimate of the true 
standard deviation a. It is intuitive that the spread in the val- 
ues of such an estimator is broader when the number of dat- 
apoints (here: radial-velocity measurements) is smaller. Let us 
assume that the radial-velocity jitter follows a Gaussian distri- 
bution, with variance (Sb) at any given Sb value. An estima- 
tor of a^, the sample variance s^jv-i (— ctq^), for an under- 
lying normal distribution follows a scaled distribution with 

— 1 degrees of freedom, where N is the number of mea- 
surements used to compute s^Ar_i. The mean of this distri- 
bution is and its variance 2a'^/{N — 1). The sample stan- 
dard deviation, sat-i (= (Tq), which is an estimator of a, in 
turn follows a scaled x distribution with — 1 degrees of free- 
dom, an expected value of c^a, and a variance of (1 — C4^)(T^. 
Here,C4 = \/2/{N - l)r(iV/2)/r((iV- l)/2) and T{x) is the 
gamma function. Therefore the standard deviation of the variable 
sn-1 for iV = 2 (as in sample I) is w 1.55 times more than for 
= 4 (a value prevailing for samples II, III, and IV). However, 
the effect of this for the upward spread of s jv-i is partially offset 
by the expected value of sjv_i being lower for N = 2 compared 
to iV = 4, by a factor of \/3/2. 

The other contributing factor is that sample I contains about 
3 times more objects than sample II. The numbers of (appar- 
ently) single objects differ similarly (728 vs. 227 using results 
from this section -see below- and from Sect. I2.4l l. which makes 
the low-probability, large-deviation tail of the ctq distribution 
more populated in sample I. At some value SAr_i.crit, the prob- 
ability of finding a value of s^-i > sjv_i^crit for any given ob- 
ject becomes less than l/(sample size). Finding such high s^^i 
is intuitively 'unlikely' - in fact, for large samples this 1 / (sample 

' Note that the prese nt figure differs siightly from the one pubiished 
bv lFamaev et aD (l2005h . which is in error. 



size) limit is equivalent to a probability of 1/e w 0.37 that no 
value of stv-i in the whole sample is found above S7v_]^ crit- 
For N — 2, this critical value, SAr_i.crit, above which the 
probability Prob(sAr_i) becomes less than 1/728 is 3.20. For 
Prob(sAr_i) < 1/227, SAr_i,crit is 2.85. The ratio of these val- 
ues is 1.12, which translates into a ratio of the ranges of s^-i 
easily populated by the two samples (purely due to their size, 
as iV = 2 was assumed for both). In fact, these two effects 
combine: an iV = 4 case with 227 objects (sample II) has 
SN-i,crit = 2.09, and comparing the ranges of sn-i populated 
by these samples gives a ratio of 3.20/2.09 — 1.53. Applying 
this factor to the dividing line in the S'fo— do diagram of sample I, 
of slope 0.57, gives an expected slope of 0.372 for sample II. 
This should be compared to an observed value of 0.3-0.4, which 
can be derived for samples II, III, IV when a dividing line going 
through point (0,0) is adopted, just as for sample I (Fig.[T^). The 
final detailed position of samples' II, III, IV dividing line is de- 
fined in a slightly different way than for sample I, namely from 
an Sb — (To diagram with the orbital motions of known binaries 
subtracted (Fig. [TJl). This adopted line has a positive intercept, 
resulting in less of a slope: 0.23S'fe+G.2. Nevertheless, the agree- 
ment between the above reasoning and the relative positions of 
the observational dividing lines is evident. 

In summary, stars falling above the dividing line may be 
safely flagged as binaries as long as 5*6 < 5 km s^^ (i.e. for 
non-Mira stars). For higher Sb values, the dividing line loses its 
diagnosis power, since Sb > 5 km s^^ is the realm of Mira and 
semi-regular variables which may have intrinsic radial-velocity 
variations with standard deviations of the order of 10 km s^^ 
(Fig. [It). For those stars, it is almost impossible to use radial ve- 
locities to detect binaries. Methods based on the proper motion 
may be used instead (see Sect. l2.2l l. 

Using the criterion based on the dividing line, 43 spectro- 
scopic binaries were identified in sample I (upper left panel of 
Fig.[T]l, corresponding to a binary frequency of 43 / 771 = 5.6%. 
The list of binary stars in sample I may be found in Table A. iQ 
of Famaey et al. (2005). 

Because the detection of binaries is difficult among Miras 
(see Sect. I2.4l l. we also compute the binary frequency for non- 
Miras (by limiting sample I to Sb < 5 km s^^), namely 38 / 
603 = 6.3%. That value will be of interest for comparing with K 
giants as discussed in Sect. l2.6l and TablelD 

2.2. Proper-motion binaries 

A few more binary stars have been found from the comparison 
of the short-term Hipparcos proper motion with the long-term 
Tycho-2 proper motion, since only the former is altered by the 
orbital motion (for orb ital periods in the range 1000 - 30000 d; 
iFrankowski et al.ll2007i) . These 'proper-motion-difference bina- 
ries' in sample I are listed in Table [T] They are not added, how- 
ever, to the spectroscopic binary frequency derived in this paper, 
because the detection efficiency of these binaries involves dif- 
ferent selection biases than those associated with the detection 
of SBs. For instance, proper-motion binaries are efficiently de- 
tected only among stars with parallaxes in excess of ~ 20 mas, 
which are rare among M giants. 



- The star HIP 26 247 = RR Cam = BP +72° 275 was erroneously 
flagged as a binary in lFamaev et all j2005h . as the result of a confusion 
between the CORAVEL coding for BD +72°275, and for the binary star 
J275 in the Hyades 
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Table 1. Supplementary 'proper- motion binaries' found in sam- 
ple I using the method described in lFrankowski et al] ( l2007h and 
extracted from their Table 2. 



HD/BD 


HIP 


Rem. 


5006 


4100 


SB" 


+58 501 


12416 


supergiant 


51802 


37391 




52554 


33929 




96572 


54571 


SB" 


103945 


58377 




108907 


60998 


SB" (P = 1714 d) 


130218 


72268 




162159 


87114 




190658 


98954 


SB" (P = 198.7 d) 






+ visual companion at 2.3" 


+39 4208 


101023 


supergiant 


199871 


103528 


SB" (P = 840 d) 


209026 


108588 




219981 


115191 




220088 


115271 


SB" (P = 2518 d) 



" the star is also flagged as spectroscopic binary from the radial- velocity 
data. 

2.3. Fast rotators ideritified from the CORAVEL line-width 
parameter Sb 

Another method of finding binaries among M giants in- 
volves identifyingfast rotators. Si ngle M giants a r e not 
expected to rotate fa st (se e e.g. iDe Medeiros et al] 119961; 
[Dc Medeiros & Mayor I l999h . so that fast rotators can be as- 
cribed to spin-up processes operating in tidally interacting sys - 
tems (like RS CVn among K giants; iDe Medeiros et al.ll2002h . 
This identification is possible by using the CORAVEL parame- 
ter Sb (measuring the average width of spectral lines corrected 
for the instrumental width) to search for outliers in the relation 
between this parameter and the stellar radius, R, as shown in 
Fig. [3l In previous papers, Sb was used as a lu minosity proxy 
dJorissen et alj 119981; IVan Eck & JorissenI 1200(1 where it was 
used to distinguish extrinsic (Tc-poor) from intrinsic (Tc-rich) 
S stars). Now we have found that the tightest correlation is in 
fact obtained between Sb and the stellar radius, rather than be- 
tween Sb and luminosity, L, or effective temperature. Toff. It is 
the Sb - radius relationship that is therefore the most appropriate 
for looking for outliers corresponding to fast rotators. 

Before embarking on the search for fast rotators, it is use- 
ful to first test the relationship between Sb and stellar luminos- 
ity, since maximum-likelihood distances are available for all the 
stars of sample I (Famaey et al. 2005). A possible relation to 
temperatu re may also b e che cked using the {V — Ic ) — Jeff cali- 
bration of lBessell et all ( IT99&) with {V - Ic) from lPlatais et alj 
(I2003h . 

Figure |2] presents the Hertzsprung-Russell diagram for sam- 
ple IJ[exx^luxling visual binaries, identified by flag 4 in Table A. 1 
of iFamaev et al] 12005), where Toff has been derived from the 
above-mentioned calibration, and the luminosity derives from 
-^^boi = Mv — (V — I a) + BCir,, w ith the bolometric cor- 
rection BCi^ from'B esselletalJ (11998') and My from Tycho- 
2 Vt2, combi i ied w ith the maximum-likelihood distance from 
iFamaev et alj (l2005l) . Two sequences are readily apparent: the 
lower, most populated one involves giant stars, and the upper 
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Fig. 2. Upper panel: The HR diagram for the full sample I 
(but the visual bin aries, identified by flag 4 in Table A.l of 
IFamaev et alj|2005l see text for an explanation about how lu- 
minosity and Tes were derived). The different symbols de- 
note Sb ranges as indicated by the labels. Lower panel: Same 
as upper panel after eliminating st ars flagged as 'Y' (young, 
supergiant stars) by IFamaev et al.l ((2005). Since maximum- 
likelihoo d distances, and henc e luminosities, could not be de- 
rived by IFamaev et al.l ( l2005h for spectroscopic binaries with 
no orbit, these binaries are displayed in both panels at an ar- 
bitrary luminosity of log L/L0 — 1.5. Evolutionary tracks from 
Charbonnel et al. ( 1996) (up to the tip of the RGB; no overshoot- 
ing) for masses 0.9, 1.0, 1.25, 1.5, and 1 .7 M^ are shown as 
solid lines. Dashed lines are tracks from iSchaller et al. I (119921) 
for masses of 2.0, 2.5, 3.0, 4.0, and 5.0 M©. For low-mass stars, 
the tip of the RGB corresponds to 5*6 = 5 to 6 km s~^. A colour 
version of this plot is available in the electronic version only. 



one, containing only stars with high Sb values, i nvolves super- 
giant stars. All stars flagged as 'Y' (young) bv IFamaev et alj 
(2005) belong to the upper sequence. The comparison of the 
evolutionary tracks with the location of the M giants in the HR 
diagram of Fig. |2| reveals that they have typical masses in the 
range 1.2 - 1.7 Mq. As far as 56 is concerned, it increases 
monotonously along the sequences; however, the existence of 
the two sequences makes it possible to check and exclude the 
correlation of Sb with either Tcff or luminosity alone. Instead, Sb 
correlates well with the stellar radius (Fig.[3]l, which is a combi- 
nation of Tcft and L. 

There are two stars (HIP 8175 = HD 10696 and HIP 98954 = 
HD 190658]] with moderate radii (log i?/i?0 = 1.6 - 1.8) and 
with Sb much larger than expected at the corresponding radii. 
At larger radii, log i?/i?0 > 2.2, the relation between R and 
Sb becomes scattered - basically because the CORAVEL cross- 
correlation dip from which Sb is derived becomes shallow and 



^ HD 10696 is not displayed on the P— 56 plot of Fig. [5] because that 
star has only one radial-velocity measurement, and was thus removed 
from sample I. 
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Fig. 3. The relationship between Sb and the stellar radius, de- 
rived from the Stefan-Boltzmann law and the HR diagram of 
Fig. [2I Open circle s denote stars flagged as 'Y' (young) by 
iFamaev et alj ( l2005l) and large crosses identify binaries. The out- 
liers with large Sb in this diagram are supposedly binaries, be- 
cause tidal interactions in the binary spun up or slowed down the 
giant's rotation. Binaries with no orbits are not included in this 
figure, since they have no maximum-likelihood distances avail- 
able either. 



asymmetric - and the identification of outliers is no longer pos- 
sible. 

HD 219654 is another example of star with too large an 
Sb, found in Fig. |2]in the form of a larger square amidst small 
ones at \o^L/Lp, = 1.5 . It does not appear in Fig. [3] be- 
cause iFamaev et al.l (|2005|) did not derive maximum-likelihood 
distances for that star; therefore, its luminosity and hence ra- 
dius cannot be derived. The reason for the lack of distance 
is that the maximum -likelihood distance estimator defined by 
IFamaev et al.l (l2005b makes use of the radial velocity, and no re- 
liable centre-of-mass velocity could be defined for those stars 
that appeared to be binaries with rather large v elocity ampli- 
tude a nd no orbit available (objects with flag '0' in IFamaev et al.l 
I2005h . 

All fast-rotator candidates are listed in Table|2] Of these, HD 
190658 and HD 219654 appear to be genuine fast rotators. The 
calibr ation of Sb in terms of the rotational velocity V sin i is 
from lBenz & MavoJ (Il981h . Fast rotation is expected in binaries 
close enough for tidal interactions to spin up the giant's rotation 
(for K giants, the critic al period below which s ynchronisation 
takes place is 250 d; see lDe Medeiros et 3112002 *). One clear ex- 
ample of this behaviour is provided by HD 190658, the M-giant 
binary with the second shortest period known (199 d; Table 8 
of Pap er I), which exhibits ellipsoidal variations as well (.Samusi 
I1997I) . Its Sb of 8.5 km s^^ formally corresponds to a rotational 
velocity of V smi = 15 km s^^. Correcting for a typical (sup- 
posedly turbulent) broadening of 3.5 km s^^ for giants of this 
radius (estimated from Fig.O yields V svai = 10 km s^^. The 
radius derived from assuming synchronisation is then R sin i — 



39 Rq. With logToff = 3.565 and log(L/L0) = 2.8, the radius 
deduced from Stefan-Boltzmann law is 62.2 Rq, which implies 
an orbit inclined at i = 39°. Rotation slower than synchronous 
would entail an inclination closer to edge-on. 

The radial-velocity data clearly reveals that HD 219654 is 
a spectroscopic binary (see its location on the Sb — aQ{Vr) 
diagram in Fig. [T^ from the parameters listed in Table |2]i. An 
outlying Sb value thus appears to be a good way to detect a 
close binary (with periods up to ~ 300 d). However, there is 
an important caveat: spurious detections are possible when a 
close visual companion contaminates the giant spectrum and 
makes it composite with seemingly broadened lines. HD 10696 
(= HIP 8175) is a such case, since it is a close visual binary de- 
tected by Hipparcos {Hp = 8.54 and 9.09 with a separation of 
0.16 arcsec). It is therefore likely that in this case the large Sb 
comes from the composite nature of the spectrum rather than 
from rapid rotation of the M giant star. 



2.4. Spectroscopic binaries (samples II, III, IV) 

The binary frequency derived for sample II, considering the sup- 
plementary knowledge gained from the extensive monitoring of 
the set of stars that had a{Vr) > 1 km s^^ (corresponding 
to sample III), and from the ao{Vr) - a{Hp) diagram (where 
a {Hp) is the standard deviation of the Hipparcos Hp magni- 
tude; see Paper I), amounts to 17/254 certain binaries, or 6.7%, 
with a binomial error of 1.6%. On top of these, 3 stars are possi- 
bly binaries, or 1.2%. To detect very long -period binaries among 
sample II, 138 of the 219 stars not present in sample III (span- 
ning the right-ascension range O'* - 16.5'' and not already flagged 
as SBs) have been re-observed, long after the last measurement 
of sample II (to make up sample IV). With this additional mea- 
surement, 5 binaries (and 2 suspected binaries) have been found 



(Table 2 of Paper I), adding a fraction of 



5x(219/138)_, 



254 



=3.1% to 



the total frequency. Therefore, the estimated fraction of spectro- 
scopic binaries among M giants is 6.7H-3.1 = 9.8%, plus 2.4% of 
suspected binaries. This frequency is 1.6 times higher than the 
value obtained for sample I, and this is related to the difference 
in binarity detection limits between samples (Sect. lZTT i. 

The above frequency should still be corrected for the se- 
lection bias acting against detecting binaries with a Mira com- 
ponent. Exact accounting for binaries among Mira variables is 
difficult. We argue that it is almost impossible to find spectro- 
scopic binaries among Mira variables (and indeed very few are 
known; ' jorissenll2003b. From the period-luminosity -radius rela- 
tionship for Miras (IWoodl[r990l) . it is possible to evaluate the 
minimum admissible orbital period for a Mira of a given pulsa- 
tion period in a detached system (assuming fundamental pulsa- 
tion and masses of 1.0 and 0.75 Mq for the Mira star and its 
companion, respectively; the conclusions are not very sensitive 
to these masses, fortunately): it ranges from about 1000 d for a 
Mira pulsating with 200 d period to a bout 2200 d for a 600 d 
pulsator (see Fig. 9.4 of |Jorissenll2003.) . These orbital periods 
translate into radial-velocity semi-amplitudes smaller than about 
10 km s^^. Given that some of these stars have intrinsic radial- 
velocity jitters of the order of several km s^^ (Fig. 1 and Alvarez 
et al. 2001), the difficulty of detecting those spectroscopic bina- 
ries using radial velocity techniques immediately becomes clear 
An alternative, astrometric method based on the comparison be- 
tween long-term Tycho-2 proper motions (not altered by possi- 
ble orbital motions ) and short-term Hipparcos proper motions 
(IWielen et al.ll 19991 and Sect. I2.2l i has not yielded any positive 
detection (because most Mira stars are quite distant, and thus 
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Table 2. Binaries inferred from the high rotational velocity of the giant primary. 



HD 


HIP 


SB?" 


CO 

(kms-^) 


Sb 

(kms-^) 


V sin i 
(kms^^) 


V - I 


Tcff 
(K) 


Mv 


R 

(Rq) 


Rem 


10696 


8175 






10.5 




2.3 


3590 


0.51 


43 


not a fast rotator; Sb contam- 
inated by close (0.16") visual 


190658 


98954 


y 


8.7 


8.5 


10 


2.1 


3670 


-0.67 


63 


companion 

Ellipsoidal variable, spectro- 


219654 


114985 


y 


6.1 


7.7 


14 


1.9 


3730 






scopic binary with P — 199 d 
not listed in Table 4b of Paper I 



because it only belongs to sam- 

plel 

" This column provides the diagnostics of spectroscopic binarity based on the location of the star in the Sb — ao{Vr) diagram (Fig.[T). 
HD 219654 : SB: JD 2448 843.634: 5.22 ± 0.4 km s'^, JD 2 449 676.270: 13.84 ± 0.3 km s"^ 



have parallax es that are too small for this method to be applied 
meaningfully: iFrankowski et al.ll2007h . 

In the Sb — ao{Vr) diagram, Mira variables have Sb > 
5 km s~^, and Fig.[T]reveals that stars with Sb> 5 km s^^ repre- 
sent a minor fraction of the total sample (less than 10%, as shown 
by the solid line in Fig.|4]i. Therefore, any correction factor will 
have a limited impact, fortunately. A simple way to estimate the 
binary frequency among M giants, circumventing the difficulty 
of detecting spectroscopic binaries involving Mira variables, is 
to restrict the sample to the 225 stars having Sb < 5 km s^^, 
for which the efficiency of binary detection is good. The binary 
frequency then becomes 0.098 x 254/225 = 0.111, to which a 
fraction of 0.027 of suspected binaries should be added. 

Finally, our samples contain a few symbiotic and suspected 
symbiotic stars: 5 in sample I (EG And, T CrB, CH Cyg, V934 
Her, and 4 Dra) and 2 in sample II (CH Cyg and 4 Dra). Their 
hot companions are most likely accreting white dwarfs, so their 
inclusion in the binary statistics of M giants is disputable, as in 
this case they correspond to the second M giant stage in the sys- 
tem's evolution. Nevertheless, EG And was flagged as binary by 
the sample I binarity criterion. If the symbiotic systems are ex- 
cluded in the evaluation of the binary frequency for sample II, 
this frequency would decrease slightly to 0.106. (Only 4 Dra has 
to be excluded, since CH Cyg, with Sb = 7 km s^^, above the 
5 km s^^ threshold, is already excluded from the above esti- 
mate.) Readers involved in binary population synthesis should 
thus note that the binary frequency of 0.1 1 1 obtained above cor- 
responds to systems with an M giant primary and a main se- 
quence or white dwarf companion. There is no guarantee, how- 
ever, that the frequency of 0.106 (obtained above by excluding 
symbiotic systems) does correspond to the binary frequency for 
systems with main-sequence companions alone. It cannot be ex- 
cluded that M giants with white dwarf companions hide among 
non-symbiotic systems, especially for M giants at the base of the 
red giant branch (RGB) where mass loss is weak. This issue will 
be developed further in Sect. 3 of Paper III. 

2.5. Binaries above Sb ^ 5 l<m 

A striking feature of Fig. [T] is the nearly complete absence of 
stars above the dashed line when 56 > 5 km s^^ in sample II. 
Sample I also exhibits an (apparent?) lack of stars above the 
dashed line beyond Sb — 5 km s^^. In fact, 4 of the 5 binaries 
observed above Sb — 5 km s^^ in sample I (Fig.[T]) are VV-Cep- 
like binaries (Cowley 1969), flagged as belonging to the young 
(' Y') group bv lFamaev et al.1 (i2005h : 

- HD 42474 (WY Gem; Sb = 8.4, ctq = 12.9 km s"!): a VV 
Cep type binary, consisting of a M2 supergiant and an early 



B type main-sequence star (iLeediaervll 1 998h ; 

- HD 190658 (V1472 Aql; Sb = 8.5, ctq = 8.7 km s'^), an 
ellips oidal or eclipsing binary of type M2.5 (iLucke & Mavoii 
fT982h : 

- HD 208816 (VV Cep; Sb = 10.2, ctq = 9.3 km s"!); 

- BD H-61°08 (KN Cas; Sb = 6.9 , ctq = 13.2 km s^^): a VV- 
Cep-type binary (ICowlevlll9^ . 

The last one, HD 219654 (V349 Peg; Sb = 7.7, ao = 
6.1 km s^^), is a giant rather than a supergiant star It is nev- 
ertheless located among supergiants because it is a fast rotator, 
spun up by tidal interactions in a close binary, as discussed in 
Sect. 1231 

To evaluate whether or not this absence of binaries above 
Sb — 5 km s^^ has any statistical significance, the Hq hypoth- 
esis that the low Sb and high 5*6 regions contain the same frac- 
tion of stars above the dashed line has been tested against the 
alternate hypothesis that the frequencies are different. The ex- 
pected number of binaries in the large Sb region (xi — P3N1) 
has been computed from the fraction of binaries in the total sam- 
ple (denoted ^3) applied to the number of stars in the large 
Sb region (Ni). The coiTesponding standard deviation on xi 
(denoted axi) is computed from the expression for a hyper- 
geometric distribution, as applies to the variable xi (the num- 
ber of binaries in the large Sb region), given the total num- 
ber of stars in the sample, N^, the total number of binaries, 
P3N3, and the number of stars in the large Sb region, A^i: 
a.i = [Nipsil ~ P3)iN3 - Ni)/{N3 - The signifi- 

cance of the reduced, squared difference — [{xi — xi) / axiY' 
is computed from the distribution with one degree of free- 
dom. Table |3] lists the (? values and the corresponding values 
of the significance level a — 0.5 Y'm\)^{x^) for different 
ch oices of the di viding line. For sample I, the stars flagged 'Y' 
bv lFamaev et al.l (|2()05) have been discarded and two versions of 
the dividing line were considered: cto = 3 km s^^ with different 
Sb thresholds and ctq = 0.575*6. For sample II, the dividing line 
was (To — 0.2356 + 0.2. Table [3] reveals that the significance 
(denoted a) of the difference is at best about a few percent. 

A Kolmogorov-Smirnov test, based on the cumulative distri- 
butions displayed in Fig. |4] has been applied to the more popu- 
lated sample I and confirms that significance level. With a total 
number of stars of 606 (after removing the visual binaries and 
the 'Y' group), there are 36 stars above the dividing line (hav- 
ing a slope of 0.57 in the 56 - ai^{Vr) diagram). The maximum 
absolute difference between the cumulative distributions of the 
total number of stars and of the number of binaries, as a func- 
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Table 3. Results of a hypergeometric test checking the Hq hypothesis that the frequency p of stars above the dividing line is the 
same on both sides of the Sb threshold value. 



Case Sb 


(70 > 3 (TO < 3 


Ni pi=Xi/Ni xi axi c' OL 


Sample I > 4.5 

< 4.5 
All 

Sample I > 5.0 

< 5.0 

A 11 

All 

Sample I > 6.0 

< 6.0 

All 


xi = 5" 154 
X2 = 17 430 
X3 = 22 584 
3 104 
19 480 

TO CO /I 

22 J 84 

3 60 
19 524 

ZZ Jo'+ 


159 0.031 5.77 2.03 0.145 0.35 
447 0.038 
606 0.036 

107 0.028 3.88 1.76 0.253 0.31 

A f\c\ r\ no o 

499 0.038 
606 0.036 

63 0.048 2.29 1.41 0.257 0.31 
543 0.035 

OUO U.UjO 




(70 > 0.57Sfe (70 < 0.5756 




Sample I > 5.0 
< 5.0 
All 


2 105 
34 465 
36 570 


107 0.019 6.36 2.22 3.85 0.025 
499 0.068 
606 0.059 




£70 > 0.2356 + 0.2 (70 < 0.2356 + 0.2 




Sample II > 5.0 
< 5.0 
All 


1 28 

29 196 

30 224 


29 0.034 3.43 1.64 2.189 0.069 
225 0.129 
254 0.118 



" Subscripts 1, 2, and 3 denote quantities for stars above the dividing line, below it, and the total number, respectively. 




I ... I . 

0.2 0.4 0.6 o.a 1 
N/N,„j 



Sb 



Fig. 4. The cumulative distributions of the total number of stars 
and of the number of binaries, as a function of Sh, for sample I 
(after removing the visual binaries as well as stars flagged 'Y' by 
iFamaev et al.ll2005l) . The inset shows how the cumulative distri- 
bution of binaries varies throughout the sample. If the frequency 
of binaries were the same for all values of Sh, the solid line 
should follow the diagonal exactly. 



tion of Sh, is 0.214 (Fig.|4|i. This difference, clearly due to the 
lack of 'binaries' with Sb > 5 km s^^, corresponds to a signifi- 
cance level of 7.5% for rejecting the null hypothesis that the two 
distributions are the same. 

The statistical significance of this lack of binaries is therefore 
not very high. Yet, such a scarcity would be easy to understand. 
First, as noted in Sect. 12.31 Sb is correlated with stellar radius. 
At larger R, for given component masses, the minimum possible 
orbital period (limited by RLOF) is longer and the maximum 
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Fig. 5. The radial- velocity jitter versus stellar radius. Circles de- 
note group 'Y' stars and crosses binary stars with orbits (binaries 
with no orbit, i.e., stars flagged '0' in Famaey et al. ( 2005) are 
not included as they do not have maximum-likelihood distances 
available, and hence radius is also not available). The curves cor- 
respond to the radial-velocity dispersion that would be exhibited 
by a giant of 1.5 Mq filling its Roche lobe in an edge-on or- 
bit (with a companion mass fixed by the mass ratio as labelled 
along the curve). The theoretical radial-velocity semi-amplitude 
has been converted into a standard dispersion by dividing by \/2. 



(edge-on) amplitude of the orbital radial-velocity variations is 
smaller. Any inclination would diminish the observed orbital 
motion even more, while not affecting the intrinsic jitter Hence, 
a binary with a more extended giant would tend to hide below 
the dividing line in the Sb — ctq diagram, among the intrinsic 
velocity jitter. This point is illustrated by Fig.|5] 
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Next, one should also notice that the RGB tip occurs around 
Sb — 5 km s~^ (see Fig.|2l), thus accounting for how the number 
of M giants drops dramatically for Sb > 5 km s^^, the RGB 
stars living longer than asymptotic giant branch (AGB) stars. 
If any, a difference in the binary properties of RGB and AGB 
stars could involve the large radii reached on the upper RGB: 
(low-mass) AGB stars would thus be restricted to long-period 
(i.e., small amplitude) systems, since the shorter-period binaries 
would have gone through Roche-lobe overflow, with dramatic 
consequences (dynamical mass transfer, common envelope, and 
orbital shrinkage to end up as cataclysmic variables) when in- 
volving a deep convective envelope, as is the case for RGB stars 
(and with none of t he physical processes to avoid orbital shrink- 
age, advocated by iFrankowski & JorissenI |2007| for RLOF in- 
volving AGB stars, being operative in RGB stars). 

This interpretation is further supported by the period - ec- 
centricity diagram (e - logP diagram) presented in Figs. |6] and 
121 revealing that no M giant is found at orbital periods shorter 
than 160 d. This limiting period corresponds to a Roche radius 
of 70 Rq when masses of 1.3 AIq for the giant and 0.6 Mq 
for the companion are adopted. In fact, no M giant binary lies 
to the left of the solid line in Fig. |6l representing the locus of 
a constant periastron distance A{1 — e) — 157 Rq which cor- 
responds to a Roche radius of 70 Rq at periastron for a system 
with such masses (see Paper III). In contrast, K giants, which are 
more compact, may be found at much shorter periods. 

However, many M giants have radii that are much smaller 
than the 70 Rq threshold obtained above: Fig. [3] shows M giants 
with radii as small as 30 Rq, in agreement with the spread ob- 
served bv lvan Belle et alj d 19991) for the radii of M giants. Stars 
with such small radii could in principle be found to the left of 
the 70 Rq RLOF limit in the e - log P diagram, but they are 
not. Part of the explanation could be that the observed envelopes 
in the e-log P diagram are defined by tidal interactions more 
than by RLOF - and tidal interactions operate well before stars 
fill their Roche lobe. This is especially indicated by the K giants, 
as they contain a prominent circularised population around and 
below the minimal period of their e-log P envelope (Figs.|6]|7]). 
Among the M giant binaries, the short-period circular subpopu- 
lation is less pronounced. Still, they stay within the 70 Rq enve- 
lope, not extending down to Roche radii of ^ 30 Rq. 

If the e-log P envelope were due to the tidal interactions, 
should it not differ in shape from the periastron RLOF type en- 
velope? For comparison. Fig. |6] also displays short- and long- 
dashed lines, corresponding to the evolution of systems during 
circularisation, which leaves the angular momentum per unit re - 
duced mass [h = A{1 - e^)] constant (IZahnlll977l:lHu3ll981h . 
Using Kepler's third law, this condition becomes P^/'^(l — e^) 
= constant in the eccentricity - period diagram. The two dashed 
lines in Fig. |6] were plotted adopting the same component masses 
of 1.3 Mq an d 0.6 Mq , and w ith the usual expression for the 
Roche radius (lEggletonlll983l) . the Roche radii corresponding 
to the circular orbits at the base of the two lines in Fig. |6] are 
70 and 200 Rq. However, lines of this shape do not follow the 
observed e-log P envelopes as closely as the lines of constant 
periastron distance do (especially in the upper panel of Fig. |7]). 
Either the shorter period objects are left outside the circularisa- 
tion line or a large empty area is included under this supposed 
envelope. This may be explained as follows. Circularisation be- 
comes fast when the giant comes close to filling its Roche lobe 
at periastron. But since the circularisation path in the e - log P 
diagram is steeper than the periastron line (see Fig. |6]l, the sys- 
tem evolves somewhat away from the periastron line and thus 
away from the strong circularisation regime. It does not "slide" 



all the way down along the circularisation path, because it loses 
the driver Only when other effects bring the star and its Roche 
lobe closer again does circularisation also accelerate again. It 
should only be fast close to the the periastron envelope. Of these 
two lines, it is thus the periastron envelope that would regulate 
the evolution of the e - log P diagram of an ensemble of sys- 
tems. 

In any case, the minimal Roche lobe radii inferred from the 
M giant e - log P diagram are apparently much too large for the 
measured radii for M giants. They correspond to Roche lobe- 
filling factors R/Rfj, which are much too small for circularisa- 
tion to operate efficiently (i.e., on a time scale comparable to the 
RGB evolutionary time scale). I nterestingly enough, a sim ilar 
difficulty has been pointed out bv lMiirset & SchmidI ( Il999h and 
iMikolaiewskal (1200 7) in the context of symbiotic stars. The for- 
mer authors have noticed that, surprisingly, the M components 
in symbiotic stars rarely had Roche lobe-filling factors over 0.6, 
whereas the latter author notes that symbiotic stars exhibit ellip- 
soidal variability despite their moderate Roche lobe-filling fac- 
tors. A possible explanation for this discrepancy could be that the 
usual expression for the Roche radius (Eggleton 1983) overesti- 
mates it due to neglecting radiation pressure and other factors 
countering stellar gravity, which can both shrink the effective 
Roche-lobe and make surface el lipsoidal distortions importan t 
at smaller radii (ISchuermanll972i: 'Frankowski & Tvlenda"200ll). 
This possibility is explored further by Dermine et al. (2009). The 
observed limit on the giant's Roche-lobe would be explained 
by a reduction of the effective gr avity on the stellar surface to 
0.15-0.35 of it s Newtonian value (IFrankowski & Tvlendall2001l: 
iDermine et al.1 '2009). 

From the above discussion one may conclude that the fre- 
quency of binaries must be expected to be lower among M gi- 
ants than among K giants, because the larger radii of the former 
forbid them from being located to the left of the dashed line in 
the e - log P diagram of Fig. |6] whereas this restriction does not 
apply to K giants. This statement will be quantified in Sect. 12.61 

2.6. Comparison with K giants 

Not many estimates of the binary freq uency among K giants ex- 
ist in the literature. An early study bv lHarris & McClurd d 19831) 
of 40 K giants in the field resulted in a frequency of 15 to 
20% spectroscopic binaries (over a 3 yr time span and with a 
radial-velocity internal error of 0.40 km s^^). In open clusters, 
Mermilliod & Mayor (2008, in prep.) find a much higher fre- 
quency of 30.8 ± 1.7% (= 217/704) for G and K giants under 
similar observing conditions (similar time span with 3 measure- 
ments per star). Surprisingly, for the giants that are not clus- 
ter members, the same authors derive a much lower binary fre- 
quency, namely 16 ± 2% (= 46/285). 

A binary frequency may also be derived for the sample of 
Hipparcos K giants monitored with CORAVEL and studied by 
iFamaev et al.l (12005 ). Table |4] lists the binary frequency among 
that sample of K giants. To compare binary frequencies among 
K and M giants, one should try to avoid systematic and selec- 
tion effects as much as possible. Most importantly, the compar- 
ison should involve binary frequencies derived from the same 
number of measurements covering approximately the same time 
span. For M giants. Table |4]reveals that not respecting this con- 
dition may cause the binary frequency estimate to vary by about 
a factor of two (from 6.3% in sample I to 11.1% in samples II, 
III, IV, after excluding the Miras). Therefore, the binary frequen- 
cies among sample I of M giants and among K giants were re- 
assessed based on the radial-velocity standard deviations com- 
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Table 4. Summary of the binary frequency among different samples of K and M giants. 



Stellar type 


Sample properties 




iVtot 


binary frequency (%) 


Ref. 


Field Mill 


I (whatever TV > 2, all Sb) 




771 


5.6 


Sect.O 




II, III, IV (N > 4, all Sb) 




254 


9.8 (+2.4?) 


Sect. [Ml 


Field Mill (without Miras) 


I (whatever > 2, S& < 5) 




603 


6.3 


Sect.O 




II, III, IV (iV > 4, S-b < 5) 




225 


11.1 (+2.7?) 


Sect. [241 


Field Kill 


(whatever N > 2) 




5643" 


14.5'' 


Famaev et al. (2005) 




Mermilliod & Mayor (2007) (TV 


= 3) 


285 


16 ±2 




Cluster Kill 


Mermilliod & Mayor (2007) (N 


= 3) 


704 


30.8 ± 1.7 





" This number differs from the one quoted by Famaev et al.. (.2005,) , because we restrict the sample to stars having at least two radial-velocity 
measurements. 

' at the level Prob(x^(V^r)) < 0.001 




p (d) 




1000 



p (d) 



Fig. 6. Top panel: The distribution in the eccentricity - pe- 
riod diagram of the 164 binaries with a Kill primary ex- 
tracted from the Ninth Catalogue of Spectroscopic Binary Orbits 
jPourbaix et al.ll2004l) . The dashed lines coiTespond to the evo- 
lution of the systems during circularisation. The thick solid line 
represents the locus of constant periastron distance A{1 — e) = 
157 i?0, corresponding to a Roche radius at periastron of 70 R©, 
assuming masses of 1.3 and 0.6 A/q for the giant and its com- 
panion, respectively. Most of the binaries with P < 100 d are 
fast-rotating, active RS CVn stars. Some of these systems are 
not yet tidally circularised. They could actually be pre-main- 
sequence stars (PMS ), since RS CV n and PMS sta rs are some- 
times confused (e.g., Rucinskilll977[) . Star 28 from [Torres et alj 
(I2002h is one such case (P = 40.9 d, e = 0.40, spectral type 
K3III). Bottom panel: Same as top for the M giants. The short- 
dashed line coiTesponds to evolution of longer-period systems 
during circularisation (see Sect. l2.5l for discussion of the lines). 



puted only from 2 datapoints (the first and the last available). 
The stars flagged as binaries in a given sample may then change. 
However, it is found that the total «Mm/7er of binaries remains the 
same (to within 2 units), both for K and M giants. Therefore, we 
may conclude that the fraction of spectroscopic binaries among 
M giants (after correcting for the difficulty of finding binaries 



Fig. 7. Top panel: The distribution in the eccentricity - period 
diagram of K giants i n clusters with a turnoff mass higher than 
1.75 Mq. Data from iMermilUod et al.1 (l2007h . The solid lines 
represent the loci of constant periastron distances A{l — e) = 79 
and 157 Rq, coiTesponding to a Roche radius at periastron of 35 
and 70 R©, respectively. Bottom panel: Same as top for the M 
giants (open circles) and K giants (filled squares) in clusters with 
a turnoff mass lower than or equal to 1.75 Mq. 



among Mira variables) is less than among K giants by a factor 
6.3/14.5=0.43. 

What could be the origin of such a difference in the binary 
frequencies among K and M giants? Part of it probably comes 
from the greater difficulty of finding binaries among M giants 
because of their smaller orbital velocity amplitudes and larger 
intrinsic jitter (Fig.|5]l, which prevents a small velocity difference 
from being ascribed to the orbital motion as it is for K giants. 
This is especially true for a small number of measurements, as 
in our sample I. 

Another possibility may reside in the different distributions 
of K and M giants in the eccentricity - period diagram, as shown 
in Fig. |6] It must be stressed that this figure refers to samples 
of field giants, i.e., orbits for K giants being ret rieved from the 
Ninth Catalogue of Spectroscopic Binary Orbits jPourbaix et al.l 
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|2004|) . ex plicitely excluding the o rbits for K giants in clusters 
derived bv'M ermilliod et al] (1200 7^ (Fig.]?]). This allows for a di- 
rect comparison of this sample with the sample of M giants with- 
out worrying about differences between field and cluster popu- 
lations. 

For stars evolving on the first-ascent giant branch, the radii 
of K giants are on average smaller than those of M giants, so that 
the minimum orbital period observed for M giants is expected to 
be somewhat longer than for K giants, as confirmed by Fig. |6] 
The resulting decrease in binary frequency may be estimated by 
computing the ratio between the number of K giants to the right 
of the solid line to the total number of K giants, or 118/164 = 
0.72. Applying this factor to the frequency of M giant binaries 
thus yields 6.3 / 0.72 = 8.8 %, closer to the binary frequency 
among field K giants, but still well below it, and this result is ro- 
bust. For example, if the long-dashed line from Fig. |6]is used as 
the limit, the correction factor would be 126/164 = 0.77 instead 
of 0.72. 

Still another bias, which may alter the factor 0.72 used 
above, should be taken into consideration. This bias is related 
to the evolutionary status of the K giant, which may either be the 
first-ascent giant branch or the core-He burning. For low-mass 
stars, the latter phase implies that stars have gone through the 
tip of the RGB where they reached a very large radius (similar 
to that of M giants). Therefore, core-He-burning K giants are 
expected to have an eccentricity - period distribution similar to 
that of M giants. 

Thus the eccentricity - period diagram of low-mass K gi- 
ants must be expected to mix such K giants with orbital proper- 
ties similar to M giants with K giants having orbital parameters 
more typical of less evolved (first-ascent giant branch) stars (i.e., 
short periods and possibly high eccentricities). Since stars spend 
more time in the He-burning clump, the former should be more 
numerous, however. For intermediate-mass K giants, this segre- 
gation does not occur since they never went through a stage with 
large radii. This may be che cked by using the samp le of K gi- 
ant binaries in open clusters (iMermilliod et al.ll2007l excluding 
stars flagged as non-members), where the K-giant mass may be 
identified with the cluster turnoff mass. Figure |2]reveals that, as 
expected, the K giants in clusters with turnoff masses lower than 
1 .75 Mq are mostly found in the region occupied by M giants, 
with only two short-period circular orbits (out of 9), in accor- 
dance with the masses of M giants estimated to fall in the range 
1.2 - 1.7 Mq from Fig.|2] The situation looks quite different in- 
deed when considering intermediate-mass K giants (upper panel 
of Fig. |7J. If the sample of K giant binaries were only made 
out of low-mass stars, the correcting factor would then be 7/9 = 
0.78 instead of 83/126 = 0.66 for the intermediate-mass K gi- 
ants (close to the value 0.72 derived above for the total sample 
of field K giants). The true correction factor should thus be the 
average of these two values, weighted by the (unknown) frac- 
tion of low-mass with respect to intermediate-mass stars among 
the sample of field K-giant binaries. This true correction factor 
should be in the range 0.6-0.8 (irrespective of the exact weigh- 
ing function), which is still too much to account for the factor 
of 2-3 difference between the binary frequencies among K and 
M giants. In the end, it is the difficulty of finding spectroscopic 
binaries among M giants due to their larger radial-velocity jitter 
that must be invoked to account for the binaries missing among 
M giants. 



3. Conclusions 

In this paper we have derived the frequency of spectroscopic bi- 
naries among field M giants, for the first time based on an ex- 
tensive sample, 771 M giants in total. The frequency obtained, 
6.3%, is much lower than for K giants (ranging from 14 to 31%, 
depending on the samples considered). However, the binary fre- 
quency derived from this large sample of M giants is not very 
meaningful, since there are important observational difficulties 
with detecting M giant binaries, because on average they have 
smaller orbital velocity amplitudes (longer periods) and larger 
intrinsic velocity jitter than K giants. This effect is especially 
important for samples with only a fewmeasurements per star, as 
is the case with this sample of 771 M giants. 

A higher binary frequency was obtained in a smaller M gi- 
ant sample with more radial-velocity measurements per object; 
11.1% confirmed plus 2.7% possible binaries. Even though this 
frequency is close to the lower bound of the binary frequency 
among K giants, we have shown that they cannot be directly 
compared because they were obtained under different observing 
conditions (mostly the number of observations per object). When 
one tries to compare the binary frequencies for samples of K 
and M giants under similar observing conditions, the binary fre- 
quency among M giants remains lower than that among K giants 
by a factor of about 2. Part of this difference may stem from the 
destruction of shorter-period K giant binaries due to catastrophic 
RLOF leading to common envelope, before they can become M 
giant binaries. However, how important this effect is depends on 
the mass distribution in the considered population. 

The lack of M-giant binaries with orbital periods shorter than 
160 d, as compared to K-giant binaries that have orbital periods 
as short as 4 d, has been clearly demonstrated from our extensive 
set of orbital elements. 

We have found that the CORAVEL line-width parameter Sb 
is better correlated with the stellar radius than with either lumi- 
nosity or effective temperature separately. This allowed identifi- 
cation of the R-Sb relation outHers HD 190658 and HD 219654 
as fast rotators, possibly due to the binary companion influence; 
indeed, HD 190658 and HD 219654 turn out to be known as 
spectroscopic binaries. A third candidate was rejected, because 
its large Sb was caused by light contamination from a close vi- 
sual companion. 

Finally, we have shown that M giants are circularised even 
though their Roche lobe-filling factor, estimated from the usual 
expression for the Roche radius, is about 0.5, in contrast to pre- 
dictions from tidal circularisation theory. A similar difficulty 
has been reported for ellipsoidal variables among symbiotic sys- 
tems. We argue that these discrepancies call for a revision of the 
Roche-lobe concept in the presence of radiation pressure from 
the giant component. 
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Fig. 1. (a) Top left panel: The radial-velocity standard deviation ao{Vr) (corrected for the average instrumental error <e>) as a 
function of the average intrinsic line width Sb, for sample I. For non-binary stars (open squares), the maximum radial-velocity jitter 
increases with Sb. Filled circles denote binary stars, as identified from their position above the dashed line, combined with the more 
detailed information gathered from samples II, III, and IV. Diamonds refer to supergiants. (b) Top right panel: Same as top left panel 
but for samples II, III and IV. Symbols are as follows: Filled circles: SBs; crosses: SB or jitter; diamond: jitter The dashed line very 
clearly marks the limit between SBs (above) and stars with intrinsic radial-velocity jitter (below), (c) Bottom left panel: Same as 
(b), adding extrinsic S stars from Udry et al.. (.1998a) . Filled triangles: non-Mira S stars; star symbol: Mira S stars, (d) Bottom right 
panel: same as (c) after removing the spectroscopic binaries with no orbit and replacing the radial-velocity standard deviation for 
orbital solutions by the standard deviation of their O-C residuals, corrected for the average instrumental error (filled circles). Filled 
triangles correspond to orbits for barium and extrinsic S stars from .Udrv et al.. (■1998a.b.) . The dashed line and open squares are the 
same in panels (b), (c), and (d). 



